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Abstract 

We study the large-* behaviour of the physical evolution kernels for flavour non-singlet observables 
in deep-inelastic scattering, where x is the Bjorken variable, semi-inclusive e + e~ annihilation and 
Drell-Yan lepton-pair production. Unlike the corresponding MS-scheme coefficient functions, all 
these kernels show a single-logarithmic large-jc enhancement at all orders in 1 —x. We conjecture 
that this universal behaviour, established by Feynman-diagram calculations up to the fourth order, 
holds at all orders in the strong coupling constant oc s . The resulting predictions are presented for the 
highest ln n (l— x) contributions to the higher-order coefficient functions. In Mellin-/V space these 
predictions take the form of an exponentiation which, however, appears to be less powerful than 
the well-known soft-gluon exponentiation of the leading ( 1 —x) 1 In' 7 ( 1 —x) terms. In particular in 
deep-inelastic scattering the l/N corrections are non-negligible for all practically relevant N. 



1 Introduction 



Disregarding power corrections, hard hadron processes are described in perturbative QCD in terms 
of process-dependent short-distance coefficient functions (mass-factorized partonic cross sections) 
and universal space- and timelike parton densities including non-perturbative long-distance effects. 
The separation between the coefficient functions and the parton densities and the splitting func- 
tions governing their scale dependence is, of course, not unique beyond the leading order (LO) 
in perturbative QCD. It is usual to perform this separation in the modified [1] minimal subtrac- 
tion [2] scheme, MS, see also Ref. [3], of dimensional regularization [4], the standard framework 
for higher-order diagrammatic calculations in quantum field theory. 

While this scheme does not provide a physical definition of the parton densities, it does lead 
to a stable (order-independent) functional form of the dominant diagonal (quark-quark and gluon- 
gluon) splitting functions in the limit of large momentum fractions x [5-10]. This feature assists 
a stable evolution of the parton densities over a wide range of scales as required, e.g., for LHC 
predictions based on data from fixed-target and HERA experiments. 

The coefficient functions, on the other hand, receive double logarithmic large-x enhancements 
in the MS scheme, i.e., terms up to (1 — jc) In (1— x) occur, for all k, at the n-th order of 
the strong coupling constant a s (the offset a > 1 depends on the observable and the power k in the 
expansion in 1— x). The highest leading (k = 0) logarithms can be resummed by the soft-gluon 
exponentiation [11] which is now known to the next-to-next-to-next-to-leading logarithmic (N 3 LL) 
accuracy for inclusive deep-inelastic lepton-proton scattering (DIS), Ip — > l + X [12], Drell-Yan 
(DY) lepton-pair production and Higgs production in proton-proton collisions [13-15], and semi- 
inclusive electron-positron annihilation (SIA), e + e~ — > h+X (h = K, K, . . .) [16, 17]. On the other 
hand, recent studies of the subleading k = 1 logarithms [18, 19] have not led to similarly systematic 
predictions for higher-order coefficient functions yet. 

An alternative description of hard processes can be obtained by eliminating the parton den- 
sities, leading to physical evolution kernels (also called physical anomalous dimensions) for the 
scale dependence of observables [3], see also Ref. [20]. This is especially simple for the flavour 
non-singlet quantities dominating the large-x limits of the semi-leptonic DIS, SIA and DY pro- 
cesses mentioned above. The direct relation between different processes via the universal parton 
densities is absent in this approach, but the soft-gluon exponentiation [11] guarantees an only sin- 
gle logarithmic k = higher-order large-* enhancement [21], see also Refs. [22]. 

Using the coefficient-function results of Refs. [23-32] for the above processes, one finds that 
the corresponding non-singlet physical kernels exhibit only a single logarithmic enhancement 
for all values of k at least to the next-to-next-to-leading or next-to-next-to-next-to-leading order 
(NNLO or N 3 LO) in the expansion in a s . We are thus led to the rather obvious conjecture that this 
behaviour, already established to order in DIS, persists to all orders in oc s . The required cancel- 
lation of double logarithms in the physical kernels then implies exponentiations also of terms with 
k > 1 , yielding explicit all-order predictions for the highest logarithms in the respective quark co- 
efficient functions. In the rest of this article we derive and discuss these predictions, emphasizing 
the subleading k = 1 logarithms. Especially for this case a formal proof of the exponentiation may 
be expected in the near future from the new path-integral approach of Ref. [19]. 
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2 Physical evolution kernels for non-singlet observables 



We start by recalling the construction and fixed-order properties of the physical evolution kernels. 
We first consider the DIS structure functions (see Ref. [33] for a general overview) 

Ji=2Fi ins> 5 2 = -F 2 , ns , J 3 =F 3 V +*. (2.1) 

x 

The longitudinal structure function Fl = F 2 — 2xFi has been addressed already in Ref. [34]. Dis- 
regarding terms suppressed by powers of 1 / Q 2 , the non-singlet quantities (12.11) are given by 

7 a (x,Q 2 ) = [C a (Q 2 )®q a . ns (Q 2 )](x) = Y^aiiQ 2 ) [c a j®q a . ns (Q 2 )](x) • (2-2) 

1=0 

As usual x is the Bjorken variable, and Q 2 = —q 2 the negative squared four-momentum of the 
exchanged gauge boson. c a j represents the /-loop non-singlet coefficient function for f a with 
Ca,o(x) = 5(1— jc). The exact three-loop results c aj 3(x) for the structure functions (12.11) have been 
computed in Refs. [26-28]. Beyond this order only the Cfti 1 ^ 1 leading-n. terms are exactly known 
[35,36]. Furthermore q a ^ ns denotes the corresponding combination of the quark densities (including 
electroweak charge factors), and (g> stands for the standard Mellin convolution, given by 



[a®b](x) = £*y a (y)b(?) (2.3) 



• l dy 

y'" 

for two regular functions and Eq. (3.4) of Ref. [21] if a +- distribution is involved. The renormal- 
ization and factorization scales /j r and /u { have been set to the physical hard scale Q in Eq. (12.21) . 

The scale dependence of the running coupling of QCD, in this article normalized as 

* S 4^ ' 

is governed by 

= PW = -Ia- P , . (2.4) 

Besides the scheme-independent Po= 11/3 Ca— 2/3 [37] (with C A = N co \ ours = 3 in QCD) and 
Pi [38], also the coefficients p 2 and P3 have been computed [39,40] in the MS renormalization 
scheme adopted throughout this study. All these four coefficients are required for calculations 
including the N 3 LO quantities c a ^. Here and below ny denotes the number of effectively massless 
flavours (mass effects are not considered in this article). Finally the evolution equations for the 
quark densities in Eq. (12.21) read 

J^2^,ns(^e 2 ) = [Pa(Q 2 )®q a ,ns(Q 2 )](x) 

= £ai +1 ^,/®^,ns(e 2 )]W- (2-5) 
1=0 

As the coefficient functions c a j, the (/ + l)-loop splitting functions P a j depend only on x for the 
above choice of /j r and /j f . All three independent third-order (NNLO) non- singlet splitting func- 
tions P a 2 (x) are known from Ref. [6]. 
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The convolutions in Eqs. (12.21) and (12.51) correspond to simple products of the respective Mellin 
transforms given by 



a N = I dxx N - l a{x) (2.6) 



o 



for regular functions such as q a ns and 



a N = j^dx {x N ~ l - 1) a(x)+ (2.7) 

for +-distributions such as the leading large-* contributions to c a j (x) . Hence calculations involv- 
ing multiple convolutions as, e.g., in Eqs. ( 12.91 ) below are best carried out in /V-space where the 
coefficient functions and splitting functions are expressed in terms of harmonic sums [41]. We 
mainly use Form [42] and TForm [43] to manipulate such expressions, to transform them back 
to the x-space harmonic poly logarithms [44], see also Ref. [24], and to extract large-x coefficients 
from the results. 

It is convenient, both phenomenologically - for instance for determinations a s - and theoret- 
ically, to express the scaling violations of non-singlet observables in terms of these observables 
themselves. This explicitly eliminates any dependence on the factorization scheme and the asso- 
ciated scale /u { , and avoids the non-negligible dependence of the MS-scheme initial distributions 
for q a ns on the perturbative order. The corresponding physical evolution kernels K a can be derived 
for = Q 2 by differentiating Eq. (12.21) with respect to Q 2 by means of the respective evolution 
equations (12.41) and (12.51) for a s and q a ,ns, and then using the inverse of Eq. (12.21) to eliminate q a ^ 
from the result. This procedure yields the evolution equations [21] 

d ^ f D / \ , at \ dC a (a s ) 



dlnQ2 r a = {P.Ca.J + PC^-^OQCa.)- 1 }®^ 

1=0 

= \a s P a ,0+Y< a i +l ( P aJ-L$kCa^k))® ?a- (2.8) 

L 1=1 \ k=o ') 

Notice that in A^-space the second term in the first line simply is |3(a s ) Jin C a (a s ) / da s . Up to N 4 LO 
(terms up to / = 4 included in the sums) the expansion coefficients c a i{x) in the last line read 



C a ,l 
C a ,2 


= CaA 

= 2c a , 2 


„®2 
c aA 




Ca,3 


= 3c £ , ;3 


- 3 C a ,2 $ 




C a ,A 


= 4c Q: 4 


-4c fli3 $ 





,®3 



-2c% + 4 Ca *®c«t-cf* 



(2.9) 

with f® 2 =f<S)f etc. The above expressions for K a j>i are valid for /u r = Q, the explicit general- 
ization to fj r 7^ Q to this order can be found in Eq. (2.9) of Ref. [21]. 

The N 3 LO physical kernels for the structure functions (12.11) are not completely known at this 
point, as the four-loop splitting functions P a ,3(x) contributing to Eq. (12.81) have not been derived so 
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far beyond the small leading-/?^ contribution [45]. Already the corresponding three-loop splitting 
functions P a ^_ , however, have only a small impact at x > 10~ 3 , see Fig. 7 of Ref. [6]. Moreover, the 
dependence of the non-singlet splitting function on N and on the specific quark combination is such 
that a single four-loop moment of any of them sets the scale for the N 3 LO contributions outside the 
small-* region, cf. Fig. 1 of Ref. [6]. Such a calculation has been presented in Ref. [46], and the 
fourth-order correction is indeed found to be small. Hence a rough estimate of P a ,3(x), for instance 
via an /V-space Pade estimate, is sufficient in Eq. (12.81) for all practical non-singlet analyses. 

The expressions for the transverse, longitudinal and asymmetric fragmentation functions 

r T = F T h ns, !Fl = Fl%, Ta = Fa (2-10) 

in semi-inclusive e + e~ annihilation (see Ref. [33] for a general overview), e + e~ — > y / Z — > h +X, 
are completely analogous to those for the corresponding deep-inelastic structure functions. The 
scaling variable in Eq. (12.21) now reads x = 2pq/Q 2 where q with q 1 = Q 2 > is the momentum 
of the virtual gauge boson, and p that of the identified hadron h. The second-order non-singlet 
coefficient functions c a ^{x) for these cases have been calculated in Refs. [29], see also Ref. [9] 
where we have derived the corresponding timelike splitting functions P a 2 for the evolution of the 
non-singlet fragmentation densities q a ,ns of the hadron h. In these cases we know the three-loop 
coefficient functions c a ,3, beyond the leading large-* terms of Refs. [16, 17], only up some terms 
involving = ft 2 /6, cf. the hadronic Higgs decay rate in Ref. [10]. These incomplete results have 
not been published. Their (complete) highest In' 1 ( 1 — x) terms will be presented in the next section. 

Finally we also consider the non-singlet quark-antiquark annihilation contribution to the total 
cross section for Drell-Yan lepton-pair production, pp/pp — ► +X, 

1 J0m 

^Y = -^- (2-11) 

In a rather schematic (but for our purpose sufficient) manner this quantity can be written as 

W*,£ 2 ) = [cDY(Q 2 )®q(Q 2 )®q(Q 2 )](x) ■ (2-12) 

Here Q 2 > denotes the squared invariant mass of the lepton pair, and the scaling variable is 
given by x = Q 2 /S where S is the squared CMS energy of the proton -(anti-) proton initial state. 
As in the case of deep-inelastic scattering, q(x, Q 2 ) represents the initial-state (spacelike) quark 
densities of the proton which evolve with the splitting functions of Ref. [34]. The non-singlet 
quark-antiquark coefficient function has a perturbative expansion analogous to Eq. ( 12.21) above, 
with G in Eq. (12.1 II) chosen such that also c DY q(x) = 8(1— x). The complete expressions for the 
NNLO contribution c DY2 ( x ) have been calculated in Refs. [31, 32]. At N 3 LO only the leading 
large-* terms, (1 — x)~ l ln"(l —x) with n = 0,. . . ,5, are presently known from Ref. [13], see also 
Ref. [14]. 

The derivation of the physical evolution kernel for jF DY proceeds completely analogous to the 
paragraph of Eq. (12.81) . with the non-singlet quark-quark splitting function occurring twice instead 
of once. As we will show in the next section, this modification is irrelevant for the purpose of 
this article, the determination of subleading large-* / large-/V double logarithms in the higher-order 
coefficient functions for the quantities in Eqs. (12.11) . (12.101) and (12.1 II) . 
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3 Known large-* logarithms at the second and third order 



We next need to address the expansions in powers of ln(l— x) of the known non-singlet splitting 
functions [6,9,47,48] and coefficient functions [23-32] in the MS scheme. The spacelike splitting 
functions (12.51) for all three types of non-singlet combinations of quark densities, 

are given by 

= + B l+l d(\-x) + C I+l \n(\-x) + o((\-x) k ^ l \n'(\-x)) . (3.2) 

(1-*)+ V J 

The last term indicates that the (/ + l)-loop splitting functions P a ,i(x) receive contributions from 
terms no higher than In (1— x), and that these contributions occur at all orders in (1—x) from the 
first. It is interesting to note that the new colour structure d abc d a b c entering the valence splitting 
function P v at three loops contributes only non-leading terms (1— x)~ In (1—x), in striking con- 
trast to its dominance in the small-jc limit [6]. As indicated the first three terms in Eq. (13.21) are the 
same for all three splitting functions P a - with the non-vanishing (Z > 0) coefficients C/+i being 
combinations of lower-order cusp anomalous dimensions A„</ - and their functional form is inde- 
pendent on the perturbative order /. This independence is established to all orders in Ref. [5] for 
the first two terms, and strongly suggested for the third term by the conjecture of Ref. [8] and its 
third-order verification in Refs. [9, 10], see also Ref. [49]. 

Eq. (13.21) also holds for the corresponding timelike splitting functions [9, 48] governing the 
evolution of the non-singlet fragmentation densities, with the same large-x coefficients except for 
a sign change of C relative to the spacelike case [8]. Hence it appears that none of the non-singlet 
splitting functions exhibits any large-* double logarithms at any order of ( 1 — x) . 

The known coefficient functions for the deep-inelastic structure functions (12.11) and the frag- 
mentation functions (12.101) - with the obvious exception of J L - receive the same highest double 
logarithmic contributions, 

cafe) = jj^y^C F ) 1 p qq {x)\n 2l -\\-x) + o({\-x) k ^- l \n 2l - 2 {l-x)) (3.3) 

with 

P qq (x) = j^j--2+(l-x) (3.4) 

and Cf = (2N c )~ l (N 2 — 1) = 4/3 in QCD. Eq. (13.31) conforms to the general observation, going 
back to Ref. [50], that the coefficient of the highest (1—x) logarithm is the negative of that of 
the highest +-distribution. Actually, this pattern also applies to the C l F ~ l {Ca, n^} ln 2/ ~ 2 (l — x) and 

C l f 2 \C\, Cpjip n 2 } ln 2/ ~ 3 (l— x) terms to / = 3, see Refs. [27,28]. Analogous results, e.g., 

Cdy,/W = j^{*C F ) l Pqq {x)\n 2l -\\-x) + o((l-*)*- _1 ln 2 '- 2 (l-*)) (3-5) 
hold for the Drell-Yan cross section (12.1 II) . 
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For the convenience of the reader, we now proceed to provide the full ln"( 1 — x) contributions to 
all known coefficient functions as far as they are relevant for our main predictions in later sections. 
The coefficient functions (12.21) for J \ read 

c l,iM = In(l-x) 2C F p qq (x) 

+ C F { p qq (x) ( -3/2 - 2 H ) - 8(1 -x) (9 + 4£ 2 ) + 3 + 3/2 (1 -x) } , (3.6) 

c lj2 (x) = (\n\\-x)AC 2 F -\n 2 {\-x)C F %)p qq {x) 

+ ln 2 (l-x) [c|{/7 M (x)(-9-14Ho) + 6 + 4Ho-(l-x)(l+2H )}] 

+ ln(l -x) [Cj \p qq {x) (-27/2 - 4H 1;0 + 24H , + 12H - I6C2) + 8H lj0 

-ll-16Ho,o-32Ho + (l-x)(-53/2-4H 1 ,o + 8H ,o + 12Ho + 8C2)} 
+ C F fr>{ Pqq W (29/6 + 4 H ) - 3 + 5/2 ( 1 -x) } 

+ C F C A {p w (x)(8/3+4Hi i o+4H 0) o-4C2)+2+(l-x)(14-4C2)} 

+ G>(C A - 2C F ) p M (-x) (8H_i, - 4H ,o) 
+ o(ln°(l-x)) , (3.7) 
c 1j3 (jc) = (ln 5 (l-x)4C^-ln 4 (l-Jc) 10/3C|(3 + ln 3 (l-Jc)2/3C F p^ w (jc) 
+ In 4 ( 1 -x) [Cj { p qq (x) ( - 1 5 - 24 H ) + 6 + 8 H - ( 1 -x) ( 5 + 4 H ) } ] 
+ ln 3 (l-x) [c|{p M (x)(-18-8H 1) o + 296/3H 0) o+54Ho-48C2)+32H 1) o 
-22 - 64H ,o - 84H - (1 -x) (33 + 16Hi )0 - 32H , - 54H - 16 ^2) } 

+ C 2 F p { p qq (x) (70/3 + 1 64/9 H ) - 8 - 4 H + ( 1 -x) ( 8 + 2 H ) } 

+ CfC A { Pqq (x) (32/3 + 8 Hi i0 + 8 H , - 1 6 £ 2 ) + 4 + ( 1 -x) (28 - 8 £ 2 ) } 

+ C}{C A -2C F ) p qq {-x) (16H_i, - 8H ,o) 
+ 0(ln 2 (l-jc)) . (3.8) 

Here and below we suppress the argument x of the harmonic polylogarithms [44] for brevity. 
Furthermore we use a slightly non-standard set of basis functions, i.e., Eqs. (13.71) and (13.81) include 

Hi,o(x) = Hi j0 (x) + ^ 2 = -mxln(l-x)-Li 2 (x) + £2 , 
H-i ;0 (x) = H_i !0 (x) + C 2 /2 = lnxln(l+x) + Li 2 (-x) + C 2 /2 

besides Ho(x) = Inx and Ho,o(x) = 1/2 ln 2 x All (modified) H-functions employed in our equa- 
tions have a Taylor expansion at x = 1, starting at order (l—x) or higher, with rational coefficients. 
Thus also all terms with the Riemann ^-function can be read off directly from our expansions. 
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The corresponding results for the structure function ?2 m (12.11) are given by 

c 2 ,iW = c 1A (x)+4xC F (3.9) 

c 22 (x) = c 1)2 (x)+ln 2 (l-x)8xCj 

+ ln(l-x) [Cj {8-x(4 + 24H )}-4xC F po-C F (C A -2C F )16x(l-^ 2 )] 

+ C} |l2-8H - (130/3 + 8Hi, -16H ,o + I6C2)*} 
-C F p {4- (50/3 + 8 H ) x] 

+ C F (C A -2C F ) |-32/(5x 2 ) (H-i,o- &/2) -32/(5x) (1 -Ho) 
+ 8/5 - 32H_ 1)0 - 16/5 H + 16^ 2 +x (236/15 + 16 [2H_i _i,o 

- H-1,0,0 + Hi,o,o - H_ 1; o + H ,o] + 104/5H - 8 ^2 - 24 £ 3 ) 
-48/5x 2 (l+Ho)+48/5x 3 (H_ 1; o-H ,o + ^/2)} , (3.10) 

c iA x ) = c l3 (x)+\n 4 (l-x)SxC^ 

+ ln 3 (1-x) [c|{16-x(8 + 48H )}-32/3jcC|po-C^(C A -2C F ) 32* (1 - &)] 

+ ln 2 (l-x) Cj |l6-48Ho-(166/3 + 32Hi )0 - 160H 0) o-24Ho + 96^ 2 )jc} 

- Cj p {20 - (158/3 + 52 H ) x} + C F $ 4x 

+ Cj{C A -2C F ) { -64/(5x 2 ) (H_i,o - - 64/(5*) (1 - Ho) 
_ 144/5 - 64H_i j0 - 32/5H + 64^ 2 +x(872/15 + 32 [2H_i _i j0 
- H_i !a o + Hi A o - H- i , + Ho,o] + 688/5 H - 48 & - 96 & H - 48 £ 3 ) 
-96/5x 2 (l+H )+96/5x 3 (H_i,o-H ,o + W2)} 

+ C F (Q-2C F )po24x(l-^ 2 )+C F (C A -2C F ) 2 32x(^-C3) 
+ 0(ln(l-x)) . (3.11) 

The differences c 2 , — c l , are, of course, the coefficient functions for the longitudinal structure 
function in DIS. Hence we have included one more order in ln(l— x) in Eqs. (13.101) and (13.1 II) . 
This order additionally includes a combination of three weight-three harmonic polylogarithms, 

H-i,-i,o(x) = H_i,_i,o(x)+H_i(xK 2 /2-C3/8, 
H-i,o,oW = H-i^o W- 3^ 3 /4 , 
Hi,o,oW = Hi, ,o(x)-^3 

besides the unmodified Ho.o.oM = 1/6 ln 3 x. The reader is referred to Ref. [24] for expressions of 
these functions in terms of the standard polylogarithms Li 2 (x) and lAj,{x). 
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To the same accuracy as Eqs. (13.61) - (13.81) for J\, the coefficient functions for can be 
written as 



£3,2$ 



c 1)1 (jc)-2C f (1-jc) 
c lj2 (x)-ln 2 (l-x) 4C|(l-x) 



(3.12) 



+ ln(l-x) [C F ! {-32H + (l-x)(18 + 28Ho-16^ 2 )} 

+ 2C F (3o(l-x)+C F C A {16Ho-8(l-x)(l+H -^)} 
-C F (C A -2C F )8p w (-x) (2H_ l5 o-H 0; o) 



+ o(ln°(l-x)) , 

= c 1;3 (x)-ln 4 (l-x)4c3(l-x) 



(3.13) 



+ In 3 ( 1 -x) [Cj { -64 H + ( 1 -x) (36 + 56 H - 32 £ 2 ) } 

+ 16/3C|(3 (l-x)+C^C A {32Ho-16(l-x)(l+H -C2)} 
- Cj (C A - 2C F ) I6p qq (-x) (2H_i, - H 0)0 ) 



+ o(ln 2 (l-x)) . 



(3.14) 



As mentioned below Eq. (13.81) . the H-functions in our expansions start at order (l—x) or higher at 
large x. Hence one can directly read off from Eqs. (13.121) - (13.141) that the coefficient functions for 
7\ and differ for x — > 1 only in terms of order ( 1 —x) . This fact was already noted in Ref . [28] . 

The third-order coefficient functions (13.81) . (13.1 II) and (13.141) receive contributions from new 
flavour classes involving the higher group invariant d abc d a bc [26-28]. The highest d abc d a bc terms 
behave as (l—x) In (l—x) for f \ and 7%, an d as (l—x) ln 2 (l— x) for J3. Their leading contribu- 
tions for the longitudinal structure function is of order (l—x) In (l—x). These terms will not be 
relevant on the level of our present analysis. The same holds for the new three-loop functions g;(x) 
which also show only a single-logarithmic behaviour for x — > 1 [27,28]. 

The coefficient functions for the transverse fragmentation function fj are related to those 
for f\ in DIS by suitably defined analytic continuations. Hence we also present their In (l—x) 
expansions relative to the results for c x } {x) in Eqs. (13.61 ) - (13.81 ): 



c M (x)+C F {l2£ 2 5(l-x)+6p w (x)Ho-6 + 3(l-x)} , 
ci , 2 (x) + In 2 ( 1 -x) Cj {7Ap qq {x) H - 12 + 6 ( 1 -x) } 



(3.15) 



+ In ( 1 -x) [C 2 {p qq (x) (8 Hi, - 28 H , - 1 8 H + 24 ^ 2 ) 



+ 22 + 40 H ,o + 1 2 H - ( 1 -x) ( 1 1 + 20 H , + 1 H ) } 
+C F p {-6p w (x)H + 6-3(l-x)} 

+C F C A {p w (x)(-8H 1; o-4H ,o)-4 + 2(l-x)} 



+ 0(ln°(l-x)) , 



(3.16) 
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CrsW = c 1)3 (jc)+ln 4 (l-x)cJ{36p OT (x)Ho-12 + 6(l-x)} 



+ ln 3 (l-x 



Ci 



x) (16Hi, - 104Ho )0 -72Ho + 48^ 2 ) 
+ 44 + 1 28 H ,o + 24 H - ( 1 -x) (22 + 64 H ,o + 44 H ) } 
+C 2 p { -232/9 Pqq {x) H + 16 - 8 (1 -x) } 



+C^C A {p qq (x) (-16Hi )0 - 8H ,o) -8+4(1-*) 
+ o(ln 2 (l-x)) . 



(3.17) 



As for the spacelike case of Eqs. (13.91) - (13.1 IK also the coefficient functions for the timelike 
longitudinal structure function J F in Eq. (12.101) will be needed to one more order in ln(l— x) 
below. Their corresponding expansions are given by 



C L,3W 



2C 



F ! 



(3.18) 



= ln 2 (l-x)4Cj 

+ ln(l-x) [C F : {-2 + 4Ho + 4x}-2C F po-8C F (C A -2C F )(l-^ 2 )] 
+ C|[-41/3-12Hi ;0 -12Ho,o + 2Ho + 16C 2 -x(2-8Ho)} 
+ C F p {25/3-2H -2jc} 

+ C F (C A -2C F ) |-24/(5x 2 ) (H_i j0 -C2/2) -24/(5x) (1 - Ho) + 118/15 

- 8 [25-1,-1,0 -H- i A o -2H -i,o -Hi A o -H_i,o] - 12/5 H 
-4^ 2 -8^2 H - 12^3 +*(4/5 + 16H_i,o-16Ho,o + 8/5Ho + 8^ 2 ) 

- 16/5 x 2 (l+H ) + 16/5 jc 3 (H_i,o-H ,o + ^2/2) 
= In 4 (l-x)4Cji 

+ In 3 ( 1 -x) [Cj { -4 + 8 H + 8jc} - 1 6/3 Cj (3 - 1 6 Cj (C A -2 C F ) ( 1 - £ 2 ) ] 



(3.19) 



+ ln 2 (l-x) 



r 3 
c F 



-35/3- 16Hi,o-16H ,o-Jc(8-16Ho) 



+ C|p {79/3-10Ho-10x} 

+ C 2 (C A -2C F ) |-48/(5x 2 ) (H_i )0 -C2/2)-48/(5x) (1 - Ho) +436/15 

— 16[2H_i _i,o — H_i,o,o — 2Ho,-i,o — Hi,o,o — H_i,o] — 104/5Ho 
-24^2- 24^ +x(-72/5 + 32H_i )0 -32Ho,o + 16/5Ho + 32C2) 

- 32/5 x 2 ( 1 + H ) + 32/5 x 3 (H_!, - H ,o + fe/2) } 

+ 2C F p 2 + C F (C A -2C F )po 12(1-^2) +C F (C A -2C F ) 2 16(^-^3)" 
+ o(ln(l-x)) . (3.20) 
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The last two equations include one more modified harmonic polylogarithm, 

H -i,o(x) = Ho -!,()(*) +Ho(xK 2 /2 + 3C 3 /2. 

The final observable in Eq. (12.101) . the asymmetric fragmentation function J a, is analogous to 
the structure function ^3 in DIS. The oc s -expansion (12.21) of its coefficient function reads 



(3.21) 



(3.22) 



= c TA (x)-2C F (l-x) , 
= c r , 2 (jc)-ln 2 (l-x)4Q?(l-x) 
+ ln(l-x) [C^{-32H + (l-x)(18 + 12Ho-16^ 2 )} 

+ 2C F (3 (l-x)+C F C A {16Ho-8(l-x)(l+Ho-C 2 )} 
-C F (C A -2C F )8p qq (-x) (2H_i )0 -H ,o) 
+ o(ln°(l-x)) , 
= c T ^x)-\n A {\-x)ACj{\-x) 
+ In 3 ( 1 -x) [Cj {-64 H + ( 1 -x) (36 + 24 H - 32 £ 2 ) } 

+ 16/3C^po(l-^)+C^C A {32Ho-16(l-x)(l+Ho-C 2 )} 
- C 2 (C A -2C F ) I6p qq (-x) (2H_ 1;0 - H ,o) 
+ o(ln 2 (l-x)) . 

The relations (13.171 ), (13.201 ) and (13.231 ) for the third-order timelike coefficient functions have not 
been presented before. These results have been obtained by extending the analytic continuations 



-A,3 



to 



(3.23) 



of Ref. [9] to terms of order a s e in dimensional regularization. At the end of this section we will 
present sufficient evidence that, for these In ( 1 —x) contributions, this analytic continuation does 
not suffer from the n 2 -problem mentioned below Eq. (12.101) . 

Finally the known coefficient functions for the Drell-Yan cross section (12.1 II) are given by 
c DYj1 (jc) = hi(l-x)SC F p qq (x)+C F {-4p qq (x)Ho-d(l-x)(16-Sb)} , (3-24) 

c D Y, 2 to = (m 3 (l-x)64C|-8m 2 (l-x)C F p )p M (x) 

+ In 2 ( 1 -x) [Cj{ - 1 24 p qq (x) H + 64 H - ( 1 -x) (64 + 32 H ) } ] 

+ ln(l-x) 



Cf {P qq {x) ("128 - 8Hi )0 + 1 12H ,o - 24 H - 64£ 2 ) - 4 
+96 [Hi, - H ,o] - 160 H + (1 -x) (4 - 48Hi ;0 + 48H ,o + 168 H ) 
+ C F p {p qq (x) (40/3 + 16H ) + 16(l-x)} 
+ C F C A [p qq (x) (32/3 + 8Hi,o + 16H 0l o - 16 £ 2 ) + 4+ 16H lj0 
+32 H + ( 1 -x) (44 - 8 Hi j0 - 1 6 H 
+ o(ln°(l-x)) . 
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(3.25) 



At the third order only the +-distributions contributions [(1— x) ln n (l— x)] + , n = 0, ...,5 are 
known so far, see Ref. [13]. 

Eqs. (13.61) - (13.231) can be employed to derive the In (1 — x) expansion of the physical evolution 
kernels (12.81) for the deep-inelastic structure functions (12.11) and the e + e~ fragmentation functions 
7t, Tt + 7l and f A of Eq. (12. 101 ). Recalling that K cun denotes the N"LO kernel for J a , one finds 



K a,o( x ) = 2C F p qq (x)+3C F 8(1 -x) , 

K aA (x) = m(l-x)^(x)[-2C F PoT8CjH ] + o(ln°(l-x)) , 

K a ,2(x) = ln 2 (l-x) j p w (x)[2C F p 2 ±12C^(3 Ho + 32C^H ,o] + o(ln(l-x)) , 

K a j(x) = ln 3 (l-x)p qq (x) [-2C F (3 3 T 44/3C^P 2 Ho-64C^PoHo,o + ^ 3 C;h OA o] 



From the NLO result K a j we have only written down the leading ln(l— x) terms. These contri- 
butions are the same for all six structure functions up to a sign change of the Ho terms between 
the DIS quantities (upper sign) and the fragmentation functions (lower sign). The non-Po terms in 
Eqs. (13.261) are the contributions of the MS splitting functions (13.21) . consequently the fourth-order 
coefficient ^/> 3 is unknown at this point, but irrelevant for our further considerations. 

The corresponding results for the Drell-Yan cross section (12.1 II) are given by 



Eqs. (13.261) and (13.271) represent our crucial observation: the physical kernels for all seven 
non-singlet observables display an only single-logarithmic large-x enhancement, at all powers of 
(1— x), to all orders in a s for which the corresponding coefficient functions are known. We con- 
sider it extremely unlikely that this pattern is accidental, and hence conjecture a single-logarithmic 
behaviour of these physical at (all) higher orders in a s , with the leading contribution showing the 
same independence on the specific structure function as in Eqs. (13.261) . In support of this conjec- 
ture we note that for the +-distribution parts of K a ^ n (including the Drell-Yan case), recall Eq. (13.41) , 
this single-logarithmic enhancement is established by the soft-gluon exponentiation as explained 
in the next section. Furthermore the all-order leading-n^ results of Ref. [36] prove the all-order 
generalization of Eqs. (13.261) for the Cf$q contributions to the DIS kernels K\ n and K2 n . 

The single-logarithmic enhancement of the physical kernels directly leads to predictions for 
the highest ln(l— x) terms of the higher-order coefficient functions. Considering, for example, 
the third (N 3 LO) line of Eqs. ( 12.91 ), one notes that the convolutions of c a ^\{x) and c a ^{x) lead to 
terms up to In (1— x). The vanishing of terms higher than In 3 (1— x) thus fixes the ln 4 (l— x) and 
ln 5 (l— x) terms of c a ^(x). In fact, exactly this reasoning, together with the absence of any ^2 
terms in Eqs. (13.261 ), provides the additional confirmation of the correctness of Eqs. (13.171) , (13.201 ) 
and (13.231) mentioned below the latter equation. 



+ o(ln 2 (l-x)) . 



(3.26) 



^DY,o( x ) 
K DY,2^) 



4C F p qq (x) + 6C F b(l-x) , 

\n(l-x)p qq (x) [-8C F po-16C^H ] + o(ln°(l-x)) , ( 
\n 2 (\-x) Pqq (x) [l6C f (3 2 + 56Cj(3 Ho + 64C^Ho,o] + o(ln(l-x)) . 



(3.27) 
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4 Soft-gluon exponentiation of the leading contributions 



The leading (+-distribution) large-* terms of the above coefficient functions can be expressed to 
all orders in oc s in terms of the soft-gluon exponentiation [11]. Switching to the Mellin moments 
defined in Eq. (12.71) . these contributions to Eq. (12.21) can be written as 

g (a s )cxp\\nNg l (X)+g 2 (X)+a s g 3 (X) + O(a 2 f(X))\ (4.1) 



C(N) 

up to terms which vanish for N^°°. Here we have used the standard abbreviation 

X = a s (3 L = ^ p L with L = \nN , (4.2) 
4tc 

and we have again put /j r = /j f = Q. By virtue of the first line of Eq. (12.81) - the logarithmic 
derivative in Af-space - Eq. (14.11) leads to the following expression for the resummed kernels up to 
next-to-next-to-leading logarithmic (NNLL) accuracy [21]: 



K a (N) = -T] a (A l a s +A 2 a 2 +A 3 a 3 ) lnJV - (l + a s + ^ a 2 )X 

v Po Po ' 



dgq,l 

dX 



- (aspo + ^p^-g^ - a 2 %^-(Xg a , 3 (X)) + 0(a 3 s (f(X)) (4.3) 

with r\ a = 2 for a = DY and r\ a = 1 otherwise, cf. the last paragraph of Section 2. Thus the 
leading logarithmic (LL), next-to-leading logarithmic (NLL) and NNLL large-/V contributions to 
the physical kernels are of the form (a s InN)", a s (a s \\\N) n and a 2 s (a s \\\N) n , respectively. Recalling 

(_i)« rin " 1 T 1 — x) 
f(N) = ±—L\n n N+0(ln n - 2 N) for f(x) = ^_ ' 

one notes that Eq. (14.31) implies that the single-logarithmic enhancement (13.261) and (13.271) holds to 
all order in oc s for the +-distribution contributions. 

In the next two sections we will provide analogous all-order results for the subleading N~ 1 In' 7 N 
contributions to the coefficient functions. These results, however, will be restricted to a tower- 
expanded NLL accuracy, see Ref . [51]. The exponents analogous to Eq. ( 14.11 ) will be given relative 
to the LL and NLL functions entering the soft-gluon exponent. For the deep-inelastic structure 
functions Ji,2,3 an d the fragmentation functions ?t+l,t,a these functions read 



X+(1-X)\n(l-X) = ^g lk (a s L) 



k=\ 



>fl,2 



[a s L) 



k=l k(k+l) 
YeAi-fli 



(4.4) 



+ 



AiPi 



\n(l-X)-^ X + \n(l-X) 
Po 



1 



k=l 



^ + ln(l-^) + -ln 2 (l-X) = £g 2 ,(a s L 
A 2 , Aip! 



YeAi -Bi 



k=\ 



5Ti(Jfc — 1) — 1 



(4.5) 
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with Qkj = 1 for k> j and Qkj = else, and Si(k) = Y!j=i 1 / 7 [51]. Here A i andA2 are the one- 
and two-loop +-distribution coefficients in Eq. (13.21) . given by [52] 



Ai = 4C F , A 2 = SC F K 



8C f 



67 
18 



?2 )Ca 



5 

- n f 
9 f . 



(4.6) 



Note that, besides these expansion coefficients and those of the beta function (12.41) . only one addi- 
tional coefficient, 

B\ = -3C F (4.7) 



enters the function g a i in Eq. (14.51) [11]. This pattern does persist at higher orders of the expo- 
nentiation [22,25]. Consequently the functions g n> \ are completely fixed by the first term of their 
respective expansion in a s , if the cusp anomalous dimension and the beta function are known to 
a sufficient accuracy. To a large extent the predictive power of the soft-gluon exponentiation rests 
on this fact: the calculation of the N"LO coefficient function is sufficient to also determine the 
N"LL resummation function g n +\, and thus two additional all-order towers of logarithms, see, 
e.g., Ref. [12]. As we will see below, however, this situation does not directly generalize to the 
non-leading large-x/large-iV terms addressed in the next two sections. 

The LL and NLL resummation exponents for the Drell-Yan cross section (12.1 II) are related to 
Eqs. (I4~4l) and (|4~51) as follows: 



£dy,i(^) = 2 Sdis,i( 2? » 

SdY.2^) = £dis,2( 2? with B i 



and y e -> 2y e 



(4.8) 



Here, as above, y e — 0.577216 denotes the Euler-Mascheroni constant. The absence of any NLL 
resummation coefficient additional to A\ and (3i is a low-order 'accident', non-vanishing coeffi- 
cients £>DY.n occur at NNLL [53] and all higher orders. 

For later convenience we finally recall the leading contributions g 31 , defined as for g x and g 2 in 
Eqs. (14.41) and (14.51) . to the NNLL resummation function in Eq. (14.11) . The universal coefficient for 
DIS and e + e~ annihilation, first extracted in Ref. [51] from the NNLO result of Refs. [23], reads 



>a,3l 



3155 22 r Anr or 22 
— - T C 2 -40C3-8C2Ye+ T Ye^ , 



367 \ „ „ 
Ye C F C A 



+ ( 1-12^ + 24^3 )cj 



247 4, 4 2 58 
27-3&+3Yc + yTfc)C^ 



f 



(4.9) 



The corresponding result for the Drell-Yan case [31, 53] is given by 



£dy,31 



/1616 



27 
224 

■ + ■ 



56^3-32^2 Ye + 

160 \ 

. C F n 



176 



1072 



27 ' 3 fe 9 / 



3 ,e 9 

(s ) C F n f 



Ye C F C A 



(4.10) 



Having collected all relevant fixed-order and soft-gluon resummation information, we can now 
turn to our new higher-order predictions. 
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5 Non-leading large-iV/ large-* terms in structure functions 



Keeping only the leading and subleading contributions, the large-/V behaviour of the coefficient 
functions (12.21) for the structure functions f a in Eq. (12.11) can be written as 



2n i 2n-l . i . 

c a , n (N) = £ c nk L k + - £ rfW L * + Of-jL 2 - 1 ) (5.1) 



with, as in the previous section, L = \nN. At the present accuracy the leading soft-gluon coeffi- 

(n) 

cients do not depend on the structure function, thus we have written c n k instead of c a k . C\ and C3 
are identical at the level of Eq. (15.11) as discussed below Eq. (13.141) - recall that N~ 2 \n a N corre- 
sponds to (1 — jc) ln a (l— x). Note that the second sum extends to In — 1, i.e., higher by one than 

in) 

the corresponding expansion for Fi analysed in Ref. [34]. Thus the highest coefficients d^ 2n _^ at 
each order n are identical also for C\ and C 2 . Recall that also the leading logarithms of the physical 
evolution kernels (13.261) to N 3 LO are the same for a = 1,2,3. Their l/N contributions are 



K a ,i 

K a ,2 

K a ,3 



N- l L 
N- l L 2 



-2p C F -16C|, 
-2(3 2 C F - 24(3 C 2 



N _ lJ? = - 2(3 3 C F -^P 2 C 2 . (5.2) 



As discussed at the end of Section 3, the vanishing of higher than single-enhanced logarithms 
in K a >n leads to relations between coefficient-function coefficients at different orders. For the two 
highest terms at all orders one finds 

<$Li = d al 7~^TTT s ^tKm* ^ 
a ' z " 1 a,L (n — 1)! (n— 1)! 

rW _ j r„ „ „ \ Qn3c" 2 3 f ( 2 ) \ 9 n2 ci^ 2 (i) cJV 1 



rf a,2«-2 - d n (C 23 -C 12 C n ) _ 3 + jtf fl2 -tf fli0 C 12 1 _ + rf fl0 



(n-3)! L fl <° 1Z J (n-2)! a ' u (n-l)! 

6n3£^ 6 w2 cf2~ 2 } (i) C^2 l 

= an inn-, — + /* 2 i — > + d — . (5.4) 

(n — 3)! (n — 2)! I ' (n— 1)! 

<2^ is independent of a, as noted above, hence we denote this coefficient by d\\ below. 0^ in (15.41) 
has been defined below Eq. ( 14.51 ), and the coefficients h\ 2 and /z 2 i in the last line are given by 

h n = c 23 -c n c n = ^pVi 2 , 

d n h 2 i = d$ - d$c 12 ■ (5.5) 

Here the second identity in the first line arises from the soft-gluon exponentiation (14.11 ) together 
with the LL and NLL expansions (l4~4l) and (l4~3i 
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A comparison with the tower-expansion [51] of the soft-gluon resummation reveals that also 
Eqs. (15.31) and (15.41) correspond to an exponential structure which can be written as 



C a {N)-C a 

N°L k 

^[dnL + d^]a s + [d { a\L + dfl]al + exp {Lh x {a s L) +h 2 (a s L) + . . . } , (5.6) 
where also the functions are defined in terms of a power expansion, 

h k (a s L) = £^ n (a s L) n . (5.7) 

k=\ 

Notice that d\ j in Eq. (15.61) is not identical to d^\ in Eq. (15.11) - the latter quantity receives a 
contribution from the expansion of the exponential. In this notation the third tower of logarithms 
is given by 

A") A ) Q n3 h ll 3 fi , l i2\ , 6n4^if 4 /. , , \ GnSSif 5 ,! 

<2n-3 = ^((^)T^+2H + ^Z4)l(*13 + *12*2lJ + ^^12 

with h\\ = c 12 = 2C>. 

The new coefficient h kn entering Eq. (15.81) (and its lower-logarithmic generalizations) can be 
determined iteratively from fixed-order information. The exponentiation (15.61) then ensures the 
vanishing of the third-highest (and lower) double-logarithmic contributions to the physical ker- 
nel at all orders in oc s . Consequently the conjectured single-logarithmic large-* enhancement of 
the physical kernel is equivalent to an exponentiation in Mellin space beyond the leading N°L k 
contributions. 

All coefficients entering Eqs. (15.31 ), (15.41 ) and (15.81) can be determined from present information. 
The corresponding coefficients of the exponent turn out to be the same for F\ and F%. They read 

h \k = 8lk for *= 1, 2, 3 , (5.9) 

fhi = g 2 i + ^Po + 6C F , (5.10) 

^22 = fe + ^Po + yPoQ-lSCl. (5.11) 

We conclude that the l/N leading-logarithmic function hi(a s L) for DIS is identical to its soft- 
gluon counterpart (14.41) . The function h2(a s L), on the other hand, receives additional contributions 
which, it appears, prevent direct predictions of g 23 etc from Eqs. (15.101) and (15.1 II) . This situation 
is analogous to that for Fi found in Ref. [34]. Hence also here the present predictivity of the 
exponentiation is restricted to the three highest logarithms at all higher orders in oc s . 

The prefactor functions in Eq. (15.61) required to this accuracy are given by the coefficient 

d n = 2C F , (5.12) 
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and for - and ^3 , recall the discussion below Eq. (15.11) - by 

1 3 

d ifl = yC F + 2y e C F (5.13) 
df} = -C F 2 (47 + 4C 2 -18Y e -4 Ye 2 )+C F C A (^-4^ + ^Y t 

~ C M^ + h)- (5 - 14) 
The corresponding coefficients for J 2 read 

4o = yC F +2y e C F (5.15) 
5$ = -C F 2 (ll9-28^-18Y e -4 Ye 2 )+C F C A (^P-20^ + yY t 



~ c Mil + h)- (5 - 16) 

Insertion of Eqs. (15.91 ) - (15.161) into Eqs. (15.31) . (15.41) and (15.81) provides explicit formulae for 
the coefficients of the three highest l/N logarithms in Eq. (15.11 ) at all orders in a s . For brevity, we 
here only present the fourth-order results, Mellin-inverted back to x-space. For <J \ one obtains 



1,4 



4 sn -yC F 4 Lj + {^C F 4 + ^c3(3oK 6 

z> k ,8(l-x) 3 [3 3 



'.V 



[188- 128^]C 4 + [12-48C2]C^C A + ^C^o + yC 2 (3 2 | > L 5 



.V 



+ 0(L 4 ), (5.17) 

where we have used the abbreviations 

T>„ = [(l-xy l ln n (l-x)] + and L x = ln(l-x) . 

The coefficients of C F Zj, Cp (3o^ and C f (3q in Eq. (15.171) are the negative of those of the 
corresponding +-distributions given (in terms of C F , C A and ny) in Eqs. (5.4) - (5.6) of Ref. [12]. 
Hence the general pattern noted below Eq. (13.41) is part of the present exponentiation and predicted 
to persist to higher orders. The corresponding result for J 2 can be written as 

16 

c i,a{ x ) = Ci^W + yQ^ 

+{[72-64^ 2 ]C 4 -32[l-^ 2 ]C^C A -yC|(3o}L x 5 + 0(L 4 ) . (5.18) 

This result, obtained from the subleading terms of the physical kernels of F\ and F 2 , is consistent 
with Eq. (16) (which also provides the coefficient of L 4 ) of Ref. [34], derived from the leading 
large-* physical kernel of the longitudinal structure function Fp. This agreement provides a rather 
non-trivial confirmation of our approach. 
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Although it is not fully known at present, it is instructive to consider also the fourth tower of 
logarithms. The corresponding generalization of Eq. (16) of [51] to the present case (15.11) reads 

j(») a J e «3 ^"r 3 t. (1) . 9«4^ir 4 f, tl , . 1 ,3 \ , ®mh" n 7 ,3 

+ 7~^yj (h u + h l3 h 2l +h u h 22 + hi n hl^j + {h n h n + \ h n h 2\ 

+ M (n-3)\ { h ^ + 2 h21 J + (n-4)\ [ h ^ +h ^) + 2(n-5)\ kn 



-(2) I e, !3 /i"r 3 e„ 4 /z"r 4 1 id erfr 2 



The additional second- and third-order coefficients in Eq. ( 15.191 ) are 

4$ = -Cj{^+n 2 -l2^- 2 ^- Ye+4C 2 Ye-31y e 2 -4 Ye 3 ) 
„ „ /12419 35 , 781 _ 11 2 



108 3 ' ~ 18 ' " ' 3 

4? = -C F 2 (^- + 47C 2 -60^ + yYe-28^Ye-39Ye 2 -4Ye 3 



„ „ /17579 13 „ 1333 11 

„ ,1699 2 y 107 2 ,\ „„„ 
- C M^T-3^^ + 3V (5 ' 21) 



and 



/■«'» = g31 +C F (l60-f C 2-36y e )-C f p„(^ + 2C 2 -^y e ) +Po 2 (^ + 1 l ye ) 

+ (c A -2c f ){c F (?il- 4 i ?2 ) + c A (H_5 ;2 )_p (» +?2 )} , (5 . 22) 



-(C A -2C f ){(80C F -8|3o)(l-? 2 ) + 16(C,-2C F )K3-?2)} (5.23) 

with g 31 given in Eq. ( 14.91) . Thus, in contrast to Eqs. (15.101 ) and (15.1 II) , the NNLL resummation 
functions h a ,3 are not the same for a = 1 and a = 2, and the deviation of their leading coefficient 
from g 31 involves ^-functions, including ^3 in the case of 
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The only other new coefficient entering Eq. (15.191) at order a 4 is ^23. This quantity can be 
constrained, but not completely fixed, from the rather obvious extension of Eqs. (15.21) for the 
(^-independent) leading l/N behaviour of the physical kernel to the next order, 



K 



ciA 



- 2 (3q C F — ^ DI§4 Pq Cf 



(5.24) 



The first term on the right-hand-side is fixed by the all-order leading-«y result for C a [36]. More- 
over the all-x expressions (13.261) strongly suggest that terms with a lower power of Po only con- 
tribute to K a 4(N) at higher orders in 1 /N. The consistency of Eqs. ( 15.191 ) - (15.241) then requires 



h 



23 



1 

8 



823 + 0P0 + 



=DIS 4 

8 



53 
18 



34 



poQ? + 72 . 



(5.25) 



As implied by the notation used above, also this coefficient of the NLL resummation function is 
the same for all structure functions (I2.ll ). The missing information for ^dis 4 is a next-to-leading 
large-n^ contribution to the fourth-order coefficient function. Since the leading large-Hy terms were 
derived more than ten years ago, and enormous calculational progress has been made in this time, 
an extension to the next order in should be feasible in the near future. We will comment on 



relations between the rational coefficients in Eqs. (15.101) . (15.1 II) and (15.251) below Eq. (16.211) . 
The resulting next contribution to Eq. (15.171 ) reads (recall L x = In (1 — x)) 



1,4 



Li 



-c 4 



1270 iAnA9 . 400. 
-— + 1424^2 + ^^3 



+ C F C A 



1576 1312 



V 9 3 



C2- 400^3 



+ C F PO^— 



+ C F 2 p 2 



=DIS 4 



+ ■ 



520 
IT 

277 



C2 



46 



20 

3- + T 



C 2 )+C F 2 C A po(y-40^ 2 



+c F p 



• 



(5.26) 



4 3 

As expected, the coefficient of C f Pq is the negative of the corresponding coefficient in Eqs. (5.7) 
of Ref. [12]. The presumed a-independent of ^dis 4 leads to a definite prediction for the ln 4 (l — x) 
term of the fourth-order longitudinal structure function, 



c 2,4 



Li 



1,4 



+ C£ (32^2- 160^ 3 ) -CpC A (8 + 224^ -208 ^) + \2Cj^ 



Li 



- Cj p (80 - C2) + 64C 2 Q 2 (C 2 - Ca) + ™C 2 F Qp (l - C2) ■ 



(5.27) 



This result completes the independent re-derivation of Eq. (16) in Ref. [34]. 

The vanishing of the double-logarithmic iV _1 ln 6 A^ contribution to K a ${N) fixes the final co- 
efficient in Eq. (15.191) . 



h 



14 



■OPo 3 



:i4 



(5.28) 



where the second equality refers to Eq. (14.41) . Thus, up to the presently unknown number ^dis 4 , the 
four highest \/N (ox (1 — x) ) logarithms for the structure functions (12.11) are fixed to all orders in 
a s . Moreover it appears obvious from Eqs. (15.91 ) and (15.281 ) that h\ (a s L) is identical to its soft-gluon 
counterpart gi (a s L) . 
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Instead of working out the corresponding all-order iV-space formalism at the next power(s) in 
l/N, we close this section on deep-inelastic scattering by presenting the fourth-order extension of 
Eqs. (13.61) - (13.141) . recall the last paragraph of Section 3 : 

c 1>4 (*) = (ln 7 (l-*)8/3C£-ln 6 (l-*) 14/3Cjp + ln 5 (1-x) 8/3C^p 2 ) p qq {x) 
+ ln 6 (l-x) [Cj{p qq {x] (-14-68/3 H ) + 4 + 8H - (l-x)(6 + 4H ) }] 
+ ln 5 (l-x) [c^{p OT (x)(-9-8Hi j0 +448/3Ho,o + 84Ho-64^)+48Hi ) o 

-22 - 96H ,o - 104 H - (13 + 24Hi, - 48H , - 84 H - 16 £ 2 ) } 
+ Cj p { Pqq (x) (41 + 3 16/9 Ho) - 10 - 32/3 H + ( 1 -x) (41/3 + 16/3 Ho) } 
+ C|C A {/7 w (x)(16 + 8Hi )0 + 8Ho,o-24C2) + 4+(l-x)(28-8C 2 )} 
+ Cj (C A -2C F ) Pqq (-x) (16H_ 1:0 - 8H ,o) 
+ o(ln 4 (l-x)) , (5.29) 
c 2i4 (x) = c 1>4 (jc)+ln 6 (l-jc)16/3jcC# 

+ ln 5 (l-jc) [C^{16-x(8 + 48H )}-40/3xC^po-C^(C A -2C F )32x(l-^2)] 
+ In 4 ( 1 -x) Cp- { 8 - 80H - (24 + 48 Hi |0 - 288 H , - 48 H + 1 60 £ 2 ) x} 

- Cj p {112/3 - (224/3 + 104 H ) x] +Cj $ I2x 

+ C^(C A -2C F ){-64/(5x 2 ) (H_ 1)0 -C2/2)-64/(5x)(l- H ) -304/5 

- 64H_i j0 - 32/5 H + 96 & +x (424/5 + 32 [2H_i _i, - H_i j0 ,o 
+H 1AO -H-i ; o + Ho ; o]+1168/5Ho-80C2-192C2Ho-48C 3 ) 

- 96/5x 2 (l + H ) +96/5x 3 (H_!, - H , + £ 2 /2) } 

+ C F (C A -2C F ) p 176/3x(l - C 2 ) +C 2 (C A -2C» 2 64x& 2 - £ 3 ) 
+ o(ln 3 (l-x)) , (5.30) 
c 3 , 4 « = c M (x)-ln 6 (l-*)8/3C 4 (l-x) 

+ In 5 ( 1 -x) [Cj { -64 H + ( 1 -x) (36 + 56 H - 32 £ 2 ) } 

+ 20/3 C^p ( 1-x) + CjC A {32H -16(l-x)(l+H -C 2 )} 

- (C A - 2C F ) 16p w (-x) (2H_i )0 - H ,o) 

+ o(ln 4 (l-x)) . (5.31) 

The ln 4 (l— x) contribution to c l 4 involves two unknown coefficients of K a ^(x), see Eqs. (13.261) 
and (15.241) . The corresponding terms in Eq. (15.311 ) can be predicted completely. However, as 
C3 — Ci does not correspond to an observable, we have refrained from writing them down here. 
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6 Results for fragmentation and the Drell-Yan process 



As discussed above, the subleading large-x / large-/V structure of the coefficient functions for the 
fragmentation functions (12.101) in semi-inclusive e + e~ annihilation (SIA) is completely analogous 
to that of their DIS counterparts addressed in the previous section. Consequently the notation (15.11) 
can be used for the present l/N coefficients as well. Also these contributions can be resummed in 
the form (15.61) . with the first four towers of logarithms given by Eqs. (15.31) - (15.81) and (15.191) . 

The coefficient functions C F and Ca are identical up to terms of order 1 /N 2 or ( 1 —x) , cf. Eqs. 
(13.211) - (13.231) above. The leading l/N logarithms of the physical kernels (12.81) are the same for 
all three fragmentation functions J l = !F T +F L , 1 T and f A (recall L = In AO, 



K a ,2 
K a A 



N^L 2 



2p C F + \6Cj , 
2p^C F + 24(3 C| 



2p 4 C F + ^ SIA4 p 3 C F 2 . (6.1) 



The first three lines derive from Eq. (13.261) . These results are identical to Eqs. (15.21) for the DIS 
kernels except for the different sign of the non-leading large-fly terms. The close relation between 
the SIA and DIS cases suggests ^ SIA4 = ^dis 4 f° r tne fourth-order generalization in the final line 
corresponding to Eq. (I5.24I ). 

The expansion coefficients of the LL and NLL contributions to the resummation exponential 
(15.61) . fixed by Eqs. (16.11) and the vanishing of higher than single-logarithmic contributions, read 

h lk = g lk for k= 1, ...,4, (6.2) 
h 2i = Sn + ^Po - 6C F , (6.3) 
h 2 2 = fe + ^Po-yPo^-lSC 2 , (6.4) 

h 2 3 = ^3 + ^o 3 +(§-% ± )PoC F -f ^Cj-llCj. (6.5) 

The coefficients (1631) - (1631) differ from their DIS counterparts (15.101) . (15.111) and (15.251) only by 
the signs of every second term in the expansion in powers of pV The first coefficients of the NNLL 



resummation function h a ^ (defined as g 3 in Eq. (14.11) ). on the other hand, are neither the same for 
the coefficient functions Cj,a and C/ , the SIA analogue of C% , nor do they show a close relation to 
their DIS counterparts (15.221) and (15.231) . These coefficients are 



4" = S3i-C f 2 (240-|c 2 + 36y e )-C f |5„(H? + 2C 2 + ^ Te )-|5 2 (^-^ye 

-(C A -2C ? ){c F (i-^ 2 ) + C A (f -§fe)-fc(-£ "fa)} , (6.6) 
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+ (C A -2C F ){(8C F +4po)(l-C2)-8(C A -2C F )(C 3 -^2)} • (6.7) 

Finally the required coefficients of the prefactors of the exponential, again obtained by expand- 
ing Eq. (15.61 ) in powers of 0C S and comparing to the results in Section 3, are given by 

d n = 2C F (6.8) 
for both coefficient functions, the same result as in Eq. (15.121) for the DIS case, 

45) = -yC F + 2y e C F (6.9) 

4 2 1 = -C F 2 (97-20C2 + 6Ye-4 Ye 2 )+C F C A (^-4^ + ^Y e 

- c Mrs + h) (6 - 10) 



42 = C F 2 (^--157C 2 + 12C3-^Ye + 20^Ye-29Y e 2 + 4 Yt 



and 



„ „ /9325 37 „ 47 Ay 11 2 

/989 2 , 25 2 7 \ 



rfff = -|c f +2 Tl C f (6.12) 



dfi' = -C|(109-36C 2 +6ye-4y e 2 )+C f Q(^p-12i; 2 + H Te 

-^(S+h) (6 - 13) 

<fo = C 2 (^-153fo + 36fe-^y e + 36i; 2 ye-25y e 2 + 4y 3 ) 
„ „ /6745 61 , 229 11 ,\ 

+ + + (6.14) 

Except for the coefficients with ^3 (and some obvious terms with Ye) there is no direct relation 
either between Eqs. ( l6\9l) - (16.141) and their DIS counterparts (15.121) - (15.161) . (15.201) and (15.211) . 

Inserting Eqs. (16.21) - (16.141) into Eqs. (15.31) - (15.81) and (15.191 ) we arrive at explicit predictions 
for the coefficients of the four highest l/N logarithms to all orders in a s , with the fourth logarithm 
including the unknown coefficient ^sia 4 of Eq. ( 16.51) . After Mellin inversion back to x-space the 
fourth-order result for fj (and [Fa, see above) read 
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-TA\ 



-TA 



o Cp L x ~\~ 

<D k) 8(l-x) 3 



y c f + y c f Po 



+ <M104 + 32 C 2 ] - [52 - 48 £ 2 ] C^Q ■ 



424 



16 



+ C F C A 



~Ci Po (f + fe) - 3 Ci (78 - f fe) + C?Q P „ - 40fc) 



9 

/964 
VT 
20 



+ 



112 



^2-400^) 



/290 



+C F Po 2 " ^) + C F (3 3 } L 4 + 0{Ll 



(6.15) 



where we have again used the abbreviations introduced below Eq. (15.171) . As for the correspond- 
ing +-distributions, see Ref. [17], the coefficients of C F Zj, Cp$oL%, C f |3qL| and C^PpL 4 in 
Eq. (16.151 ) are the same as in Eqs. (15.171 ) and (15.261 ) for the DIS case. 

The corresponding predictions for the total fragmentation function fi = Jp + 7l lead to the 
following results for the longitudinal fragmentation function Ji\ 

c M W = ^C;L x 6 + {[36-32^ 2 ]c4-16[l-^]CjC A -^C^(3o}L x 5 

+ {c; (64^-80^3) ~CjC A (4 + 112^2- 104^) -C^po(40-^ 2 ) 

88 

+ 6C F 2 (3 2 + 32C F 2 C A 2 (C 2 -C 3 ) + yC F 2 C A (3o(l-^)}^ + 0(L 3 ) . (6.16) 

Besides an overall factor of two arising from the different definitions of Ji in SIA and DIS, this 
expression differs from its counterparts (15.181) and (15.271) for the longitudinal structure functions 
in DIS only in the coefficient of L^Cp ln 4 (l— x). Eq. (16.161) can be derived also via the physical 
evolution kernel for the longitudinal fragmentation function, in complete analogy with the DIS 
case in Ref. [34]. In fact, Eqs. (20) - (22) of that article hold for the present case as well, with 
the above difference arising from the second-order prefactor to the resummation exponential. This 
close relation between the spacelike and timelike cases does not persist at higher orders in (1— x), 
as can be seen already by comparing Eqs. (13.91) and (13.181) . 

We now turn to the corresponding results for the non-singlet Drell-Yan cross section (12.1 II) . 
The leading l/N contributions to its physical kernel are given by G 



^DY,l 


N~ 


l L 


-8p C F - 


- 32C 2 , 


Km,2 


N- 


l L 2 


-16p 2 C F 


-112p C 2 , 


#DY,3 


N~ 


l I? 


-32(3 3 C F 


+ ^DY 3 Po 


^DY,4 


N- 


1/4 


-64p 4 C F 


+ ^DY 4 Po Cf 



(6.17) 

Here the first two lines follow from Eqs. (13.271) , while the third and the fourth are the obvious 
generalization to order a 4 and a^, respectively, exploiting the complete analogy to the DIS and 



22 



SIA cases discussed above. Also these parts of Eqs. (16.171) are of some interest despite the unknown 
subleading large-Po terms. 

This can be seen from the resulting coefficients of the LL and NLL resummation exponents, 



h 2 \ 



8lk for k = h 
821 + Po + VC F , 



h 



22 



h 



23 



822+ gPo 
#23 + Po - 



7 

3 + 24 



DY 3 

24 

DY 3 



PoC F 

^DY 4 

32 



49 



C 



DY 3 



24 



r r 2 + 343 r 3 



(6.18) 
(6.19) 

(6.20) 
(6.21) 



We note that, both here and in Eqs. (15.1 II) and (15.251) for the structure functions in DIS, the coeffi- 
cients of Cp in h 2n are given by l/n times the n-th power of the corresponding coefficient in /z 21 . 
Furthermore the coefficients of fioCp m Eq s - (15.251) and (16.211) are the products of the respective 
C F and |3oC F coefficients in h 0i and h 22 - These relations seem to point towards a general structure 
for the functions h2(a s L) in Eq. (15.61 ) which, we hope, can be uncovered in some more deductive 
approach to the \/N contributions to the coefficient functions. 

The prefactor coefficients relevant for the highest three logarithms read 



d 



4) 

DY,1 



= d\] 



8C 



F ) 



d 



(1) 

DY,0 



"DY,1 



8y e c F , 

,2> 



(6.22) 



/884 



88 



- C 2 (156 - 128 C 2 - 56y e - 64 y 2 ) +C F C A ^- — 16^2 + y Y, 
„ /176 16 



(6.23) 



Together with Eqs. (16.181) - (16.201) these results lead to the third- and fourth-order predictions 



C DY,3 ( x ) — C DY,3 



2)4,8(1-*) 

f <j [2272 + 3072^]^ 



512CJ L x 5 + { 1728Cj + — Cj $q}l a x 
544 



L 3 



512^2 



CpC A 



2944 4 
— ^ + - 



DY 3 



C F 2 p -yC F P 2 }L x 3 + 0{Ll) 



(6.24) 



and 



-DYA\ 



DY,4 



4096 



©jfe,S(l-x) 



19712 



3584 



+ ^ [9088 + 20480 £ 2 ] - [1408 - 3072 £ 2 ] c/C A 



20864 | 8^ D Y 3 



Cfa-^C}tf}l* + <9(L X 4 ) , (6.25) 



where the respective third logarithms depend on the presently unknown quantity ^dy 3 - Also in 
Eqs. (16.241) and (16.251) the coefficients of the highest +-distributions and powers of L x = ln(l— x) 
for each colour factor are equal in magnitude but opposite in sign. 
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Finally we provide the generalizations of Eqs. (13.151) - (13.251) to the next order in a s . For the 
fragmentation functions (12.101) these are given by 



'TA 



(x) = c 1)4 (jc)+ln 6 (l-x)C# {32p ot (jc)Ho-8 + 4(1-jc)} 

+ ln 5 (1-x) [c A F {p qq (x) (16Hi, - 152H ,o - 108H + 48 £ 2 ) 
+ 44 + 176H , + 24 H - ( 1 -x) (22 + 88 H 0;0 + 68 H ) } 
+Cj p { -428/9 pqq (x) Ho + 20 - 10 (1 -x) } 
+C|C A {p w (x)(-16Hi,o-8H ,o) -8 + 4(1-*)} 
+ o(ln 4 (l-x)) . 

c LA {x) = In 6 (1-x) 8/3 C A F 

+ In 5 ( 1 -x) [Cj { -4 + 8 H + Sx} - 20/3 C F p - 1 6 Cj (C A - 2 C F ) ( 1 - ^ 2 



(6.26) 



+ In 4 (1-x) 



8H lj0 -8Ho,o-4Ho-32C2-x(12-16H )} 
+ C^(3 {112/3-56/3H -56/3jc} 

+ Cj(C A -2C F ) |-48/(5x 2 ) (H-1,0 - C2/2) -48/(5*) (1 - H ) +212/5 

-16[2H_ 1: _ 1: o-H_ 1A o-2H 0) -i ) o-Hi ) o,o-H_i j o]-184/5Ho-40C2 

- 24^ + 16 £ 2 Ho +x(- 152/5 + 32H_ 1)0 -32Ho,o + 16/5 H + 48£ 2 ) 

- 32/5x 2 (l + H ) + 32/5x 3 (H_i, - H ,o + k/2) } 

+ 6C2p 2 + c2(C A -2C F )p 88/3 (l-C2)+c2(Q-2C F ) 2 32(C2-C 3 ) 
+ o(ln 3 (l-x)) . (6.27) 
CaaW = c r . 4 (x)-ln 6 (l-x)8/3C 4 (l-x) 

+ In 5 ( 1 -x) [Cj {-64 H + (1 -x) (36 + 24 H - 32 C2) } 

+20/3C|p (l-*) +C|C A {32H -16(l-x)(l+H -C 2 )} 
-C|(C A -2C F )16p M (-x)(2H_ 1)0 -H 0) o) 
+ o(ln 4 (l-x)) . (6.28) 
The corresponding result for the third-order Drell-Yan coefficient function reads 

c D y,3(*) = (ln 5 (l-x) 192C 3 -801n 4 (l-x)C 2 pV) P qq {x) 

+ ln 4 (l-x) [Cj {-648/? M (x)H + 384H -192(l-x)(2 + Ho)}] 

+ (In 3 (1-x)) . (6.29) 

Unlike Eq. (16.171 ), the fourth-order generalization of Eqs. (13.271 ) involves more than one unknown 
coefficient, hence we have not included the incomplete ln 3 (l — x) contribution in Eq. (16.291) . 
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7 Numerical illustrations 



We close by briefly illustrating the numerical size of the known and new subleading large-TV contri- 
butions to the coefficient functions. For = 4 the corresponding expansions of the two- and three- 
loop coefficient function for F 2 , the practically most important structure function, are given by 

c 2:2 (N) = 3.556L 4 + 26.28L 3 +40.76L 2 -67.13L-157.3 

+ A^ 1 (7.111L 3 + 92.76L 2 + 239.5L + 214.0) + O (N~ 2 ) , (7.1) 

c 23 (/Y) = 3.160L 6 + 44.92L 5 +238.9L 4 + 470.8L 3 -620.2L 2 - 1639L-3586 
+ /V~ 1 (9.481L 5 + 21 1.9L 4 + 1393 L 3 +4157 L 2 + 5200 L + 5230) 
+ (N~ 2 ) . (7.2) 

In Fig.Q]these approximations, with and without the 1 /N terms, are compared to the exact results 
of Refs. [23,24] and [27]. At both orders the latter contributions are relevant over the full range 
of /V shown in the figure, while terms of order 1//V 2 are sizeable only at N < 5. Note that the 
classification as TV and TV -1 terms does not reflect the numerical behaviour for the TV-values of 
the figure. E.g., the third-order increase due to the ln k N contributions in the first line of Eq. (17.11 ) 
strongly resembles a linear rise, and the sum of the 7V~ ! In N terms in the second line almost looks 
like a constant. In fact, the decrease of this contribution towards large /V is very slow: only at 
/V = 1.5 • 10 2 has it fallen to half of the value at its maximum at /V = 6.6. The situation for the 
corresponding third-order coefficient functions for F\ and F3 [28], not shown here for brevity, is 
similar except at small /V where in both cases the sum of the TV and TV -1 terms is close to the exact 
result even down to N = 1, the lowest value of N used in the figures. 

The pattern of the coefficients is rather different for both the N° and TV -1 contributions to the 
corresponding coefficient functions for the transverse fragmentation function Ft, 

c T2 (N) = 3.556L 4 + 25.69L 3 + 105.6L 2 + 104.3L 

+ /Y" 1 (7.111L 4 -29.02L 4 - 111. 4L- 504.0) + (/V~ 2 ) , (7.3) 

c Ti3 (N) = 3.160L 6 + 43.34L 5 +309.3L 4 + 1017L 3 + 2306L 2 + 2090L + 9332 

+ /V~ 1 (9.481L 5 — 10.17L 4 — 362.1 L? — 3241 L 1 ) + O (N~ l L) . (7.4) 

These expansions are shown in Fig. |2] together with the exact second-order result of Refs. [29,30]. 
As adequate for an observable measured in particular at scales not too far from the Z-mass, the 
results refer to = 5 effectively light flavours. All N° contributions are positive in Eqs. (17.31) and 
(17.41) . yielding a larger soft-gluon enhancement than in the DIS case especially due to the lower 
powers of ln/V as already discussed in Ref. [17]. On the other hand, the N~ l ln k N coefficients 
change sign here, again in contrast to Eqs. (17.11) and (17.21) . This leads to smaller and negative 
l/N corrections which do not exceed 10% except for N < 1 in the two-loop case in the left part 
of the figure, where their inclusion results in a good approximation down to /V ~ 2. At the third 
order the 7Y _1 ln/Y and TV -1 contributions are not yet known. One may expect similarly relevant 
small-TV corrections from these terms to the corresponding curve shown in the right part of the 
figure. Similar results are found for the integrated and asymmetric fragmentation functions. 
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Figure 1 : The second- and third-order non- singlet coefficient functions for the structure function 
as defined in Eq. (12.21 ) in Mellin-iV space. The leading and subleading large-Af contributions (17.11 ) 
and (17.21) are compared to the exact functions for ny = 4 light flavours. The results are multiplied 
by suitable factors compensating our small choice a s = oc s / (47t) of the expansion parameter. 
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Figure 2: As Fig.Q], but for the fragmentation function Ft at = 5. Neither the exact three-loop 
result nor the corresponding coefficients of the iV -1 In iV and ./V -1 terms are known at present. 
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We now turn to the four-loop predictions derived in the previous two sections, again focusing 
on the same two DIS and SIA observables. The known and new contributions at this order read 

c 2 4 (N) = 2.107L 8 +48,71L 7 +477. 9L 6 + 2429 L 5 +5240 L 4 - 1824 L 3 - 30308 L 2 

+ (L) + N-\ 8.428L 7 + 284.3L 6 + 3324 L 5 + [18884 + 30.86 £ DIS J L 4 ) 

+ o(N^L 3 ) (7.5) 

and 

c TA {N) = 2. 107 L 8 + 46.60 L 7 + 5 14. 1 L 6 + 3126L 5 + 1 1774L 4 + 23741 L 3 + 46637 L 2 
+ 0(L) + N~ l ( 8.428 L 7 + 32.47 L 6 - 448.1 L 5 - [7315 + 26.12^ SIA J L 4 ) 

+ o(N~ 1 L 3 ) . (7.6) 

As above, the results for F 2 in Eq. (17.51) are given for = 4, and those for Ft in Eq. (17.61) for 
rijr = 5. The N° coefficients have been presented already in Tables 1 of Refs. [12, 17]. The \n 2 N 
term in both equations includes a small contribution A4/2 from the fourth-order cusp anomalous 
dimension for which we have used the respective Pade estimates of 4310 for ny = 4 and 1550 for 
rijr = 5 [12]. The fourth N~ l logarithms receive small contributions from the presently unknown 
(and most likely identical) fourth-order coefficients ^dis 4 an d ^sia 4 of Eqs. (15.241) and (I6.ll ). Values 
expected from the latter equations contribute less than 2% to the coefficients of N~ l \n 4 N. 

The presently unknown lower-fc N~ l \n k N terms can be expected to enhance the 1 /N effects 
shown in Fig. [3] Yet already now one can conclude that the pattern of the previous two orders 
appears to persist to order a 4 , e.g., that the iV -1 contributions are small for Fj at least at N > 10. 
We stress that this figure does not intend to present the best approximation to dominant N° con- 
tributions, but simply illustrates the effect of the known terms as given in Eqs. (17.51) and (17.61) . 
Rough estimates of the missing coefficient of lniV can be obtained by expanding the soft-gluon 
exponential (14.11) or (also for the non-logarithmic N° terms) via the Mellin transform of the known 
seven +-distributions given in Eqs. (5.4) - (5.10) of Ref. [12] - note that there are some typos in 
the first archive and journal versions of this article - and in Eq. (32) of Ref. [17]. The latter article 
includes also the c 2 4 — c T 4 difference of the lnAf and coefficients in Eqs. (17.51) and (17.61) . 

Finally the large-/V expansion of the second- and third-order coefficient functions for the non- 
singlet (quark-antiquark annihilation) Drell-Yan cross section (12.1 II ) are given by 

c 2)2 (N) = 56.89L 4 +185.9L 3 + 428.6L 2 + 267.6L + 442.8 

+ AT~ 1 (113.8Z 3 + 378.4L 2 + 577.3L + 53.43) + (N~ 2 ) , (7.7) 

c 23 (N) = 202.3L 6 +1282L 5 +4676L 4 + 8172L 3 + 11404L 2 + 6395L + O(l) 

+ A^ 1 (606.8L 5 +4267L 4 + [12164-4.543^ DY3 ]L 3 ) + O (A^L 2 ) (7.8) 

for rij = 5. These expansions are shown in Fig. |4] together with the exact two-loop results of 
Refs. [31,32]. The higher-order corrections are much larger in this case than in DIS and SIA. Also 
here the \/N contributions appear to be numerically rather unimportant, a feature that appears to 
persists to even lower values of iV than for the fragmentation functions. 
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Figure 3: Large-Af contributions to the fourth-order non-singlet coefficient functions for F2 in DIS 
(left) and Ft in SIA (right). Shown are the known and Af _1 contributions as given in Eqs. (17.51) 
and Eqs. (17.61) . The results have been multiplied by 25000 ~ (47l) 4 for display purposes. 
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Figure 4: As Fig.[2l but for Drell-Yan cross section (12.1 II) . using the expansions (17.71) and (17.81) with 
^dy 3 = — 400. Besides the coefficients mentioned in the caption of that figure, also the third-order 
constant-Af and iV -1 ln 2 Af coefficient are unknown in this case. 
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8 Summary and outlook 



We have analysed the \n k (l — x) contributions to the physical evolution kernels for - including the 
results already presented in Ref. [34] - nine flavour non-singlet observables in inclusive DIS, semi- 
inclusive e + e~ annihilation (SIA) and Drell-Yan lepton-pair production. It turns out that all these 
kernels include only single-logarithmic higher-order corrections, up to cx"(l — x) ln n_1 (l — x), at 
all powers k of ( 1 — x) . On the other hand, the coefficient functions from which these kernels are 
constructed received double-logarithmic contributions up to cx"(l —x) k ln 2,1-1 (1 —x) at all orders. 
This difference implies that the terms cx"(l— x) k \n ! (1 — x) with n < I <2n are functions of lower- 
order terms, i.e., a general resummation of the double-logarithmic terms at all powers of ( 1 —x) . 

The above pattern is established to all order in a s by the soft-gluon exponentiation of the 
(1 — ln^(l — x) contributions to the coefficient functions [11-17]. All-order results underpin- 
ning it at all powers in ( 1 — x) are presently known only for the leading large-n^ contributions to DIS 
structure functions [35,36]. However, all available fixed-order results on higher-order coefficient 
functions [23-32] are consistent with the behaviour described in the previous paragraph. 

In our view it is most unlikely that this consistency is accidental, given the large number of 
observables and the depth of the perturbative expansion reached especially in DIS and SIA - for 
the latter this article includes some new third-order results. Moreover it should be noted that the 
resummation of Fl in both DIS and SIA can be consistently constructed each via two different 
physical kernels: that for these quantities themselves (starting with (1— x)~ l ) and via the differ- 
ence (starting with ( 1 — jc)°) of the respective kernels, K2 — K\ and K[—Kj, for the structure func- 
tions F\z and fragmentation functions Fjj where Fj is our notation for the total (angle-integrated) 
fragmentation function. We thus definitely expect that we are observing a genuine feature of the 
coefficient functions and expect that a more deductive approach, such as that pursued in Ref. [19], 
can provide a formal proof in the near future at least for the next power in ( 1 — x) . 

We have employed the conjectured single-logarithmic enhancement of the physical kernels to 
derive the explicit jc-dependence of the coefficients of the three highest powers of the fourth-order 
DIS and SIA coefficient functions, while in the Drell-Yan case we are restricted to two logarithms 
at order oc s 3 . For this purpose we have employed a modified basis (required far beyond the weight-3 
functions shown in the article) for the harmonic polylogarithms. An extension of these results to 
higher orders in a s is possible but not necessary at present in view of the discussion given below. 

For the subdominant (except for Fl) logarithms with prefactor ( 1 — x) we have cast our results 
in the form of an exponentiation, akin to that of the (1— x)^ 1 soft-gluon effects, in Mellin-/V space 
where these terms behave as N~ l \n k N. One more logarithm can be effectively predicted in this 
case, as the one unknown parameter turns out to be numerically suppressed. Our resummation 
of the \/N terms is, nevertheless, far less predictive than the soft-gluon exponentiation (which 
predicts seven of eight fourth-order logarithms in DIS and SIA) for two main reasons: Firstly the 
prefactor of the exponential is of first instead of zeroth order in oc s , thus one more order needs to be 
calculated in order to fix the same number of coefficients. Secondly, while the leading-logarithmic 
function, usually denoted by g^a^XnN), in the exponent is the same as in the N° soft-gluon case, 
this does not hold for the higher-logarithmic functions which have an (at least presently) not fully 
predictable power expansion (from g 2 ) and do not show any universality (from g 3 ). 
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Finally we have illustrated the numerical size of the \/N contributions. It turns out that, in 
the restricted iV-region of practical interest, the logarithms at the third and higher orders essen- 
tially compensate one power of N, i.e., the N° terms together resemble a linear increase with N, 
and l/N corrections almost look like a constant. The sum of all N° and A^ 1 contributions is 
found to provide an excellent approximation of the exact results, except at small iV- values such as 
N < 5, wherever both are known. However, only in the DIS case do the 1 /N terms constitute a 
phenomelogically significant correction over a wide range of moments. 

The main application of the present results and, hopefully, their future extensions in a more 
deductive approach - we note that also an extension of Ref. [36] to the next-to-leading large- 
rif terms would provide very useful information in the present context - may be in connection 
with future higher-order diagram calculations, e.g., of the fourth-order DIS coefficient functions: 
Firstly they can serve as important checks of such computations which will be of unprecedented 
complexity. Secondly, they will be very useful in combination with future partial results such as 
a fourth-order extension of the fixed-/V calculations of Refs. [54], as fewer computationally costly 
moments will be required for useful x-space approximations along the lines of, e.g., Ref. [21]. 
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